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Abstract 

We use chiral perturbation theory to study the extraction of pion electromagnetic polarizabilities 
from lattice QCD. Chiral extrapolation formulae are derived for partially quenched QCD, and 
quenched QCD simulations. On a torus, volume dependence of electromagnetic observables is 
complicated by SO{A) breaking, as well as photon zero-mode interactions. We determine finite 
volume corrections to the Compton scattering tensor of pions. We argue, however, that such 
results cannot be used to ascertain volume corrections to polarizabilities determined in lattice 
QCD with background field methods. Connection is lacking because momentum expansions are 
not permitted in finite volume. Our argument also applies to form factors. Volume effects for 
electromagnetic moments cannot be deduced from finite volume form factors. 
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I. INTRODUCTION 



Electromagnetic polarizabilities encode fundamental properties of bound states. The 
electric polarizability of the ground state hydrogen atom, for example, a E = aj s N /m e -Eg, 
represents the ease at which the atomic electron cloud deforms in an applied electric field. 
Here aj s = e 2 /47r is the fine structure constant, m e is the electron mass, E is the ground 
state energy, and N is a pure number, which turns out to be 9/8. Atomic polarizability data 
are well described by theoretical calculations using atomic wave-functions of the weakly 
bound electrons. Hadronic polarizabilities, on the other hand, involve non-perturbative 
physics. The electrically charged quarks inside hadrons respond to applied electromagnetic 
fields but against the strong chromo-electromagnetic forces that confine them into bound 
states. If the pion were a weakly bound system of quarks with mass m q , we might expect 
its electric polariazability to be of the form, a\ ~ aj s N /m q m^. The actual behavior is 
considerably different, a^ E = af s N/m 7T A'^ c , where A x is the chiral symmetry breaking scale. 
It thus appears that the pion cloud of the pion is what deforms in the applied field, and 
that the relevant energy scale is A x , which is an order of magnitude greater than the pion 
mass. Compared to the weakly bound scenario, the electric polarizability is a few orders of 
magnitude smaller, which indicates stiffness of quarks inside hadrons. 

Chiral peturbation theory (%PT) [l[ provides a low-energy effective theory of QCD from 
which the pion polarizabilities can be calculated in terms of a few low-energy constants 0] ■ 
At leading order in the chiral power counting, calculated values for the pure number N are 
N n = — 1/3 for the neutral pion, and N n « 1/6 for the charged pions. Comparing these 
polarizability predictions to experimental data is unlike the situation with atomic polar- 
izabilities. Without stable targets, experimental determination is considerably challenging 
at best. Pion polarizabilities, however, have been probed indirectly in several experiments. 
Three reactions are used: radiative pion-nucleon scattering (irN — > irN'y), pion photopro- 
duction in photon nucleus scattering (jA — > 7 An), and pion production seen in electron- 
positron collisions (7*7 — > nn). Neutral pion polarizabilities have been accessed only by the 
last reaction by the Crystal Ball Collaboration (3)]. The most recent experimental effort has 
been by MAMI at Mainz in measuring the difference of electric and magnetic polariz- 
abilities of the charged pion through radiative pion-nucleon scattering, and by Compass at 
CERN [|3| to measure charged pion polarizabilities using photoproduction off lead. In the 
latter experiment, final data are being taken, and soon will be analyzed. After the 12GeV 
upgrade, Jefferson Lab has plans to measure pion polarizabilities in the future. 

Experiments to determine pion polarizabilities have one feature in common: disagreement 
with predictions from chiral perturbation theory. Considerable effort has been expended 
to determine polarizabilities to two-loop order in %PT 0, 0, Isl, 0], but discrepancy with 
experiment remains. Because these experiments are indirect, the challenge is removing 
the hadronic background to isolate the signal. This is a largely model-dependent process 
with uncontrolled systematic error. Recent dispersion relation calculations, however, appear 
consistent with experimental values for the polarizabilities (lol 11]. Thus it remains unclear 
whether disagreement between theory and experiment has its roots in the experimental 
analysis, or in the behavior of the chiral expansion. 

As a first principles method, lattice QCD [l2| can be employed to determine pion po- 
larizabilities. Currently and foreseeably this is itself a considerable challenge, but progress 
has been made with background field methods [H, 0, [H, M, H H H 0, H Q . Such 



calculations suffer a number of systematic errors, such as: quenching or partial quenching, 
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quenching of sea quark charges, and volume effects. While predictions of physical polariz- 
abilities are not currently possible from lattice QCD alone, forthcoming lattice QCD data on 
polarizabilities can be used as a diagnostic for %PT. The predictions of %PT can be tested 
against lattice QCD data. To this end, we perform a one-loop analysis of the quenching and 
partial quenching effects, as well as the volume dependence of pion Compton scattering. As 
polarizabilities are the coefficients at second order in an expansion in photon momentum 
u), one would naively expect that finite volume corrections to polarizabilities can be deter- 
mined from momentum expanding the finite volume Compton tensor. We find this is not 
the case. There are many terms in the finite volume Compton tensor not anticipated by 
infinite volume gauge invariance. All terms, moreover, are form factors in uL, where L is 
the spatial size of the lattice. Because of momentum quantization, these form factors cannot 
be expanded in uL. Thus the infinite volume connection between the frequency expansion 
and the polarizabilities is lost. As polarizabilities are typically calculated in lattice QCD 
using background field methods, this means we cannot use the finite volume Compton ten- 
sor to deduce finite volume corrections to polarizabilities extracted from background field 
correlation functions. The same problem exists for electromagnetic moments. Their volume 
effects cannot be deduced from series expanding finite volume electromagnetic form factors 
about zero momentum transfer. 

Our work has the following organization. First in Section III Al we detail our conventions 
for Compton scattering and the electromagnetic polarizabilities of the pion. In Section III Bl 
we review the low-energy effective theories of QCD, and partially quenched QCD. Quenched 
QCD is discussed in Appendix |A] These theories are then utilized in Section Hi Cl to compute 
the pion electromagnetic polarizabilities in infinite volume. Next in Section lITTj we consider 
the modifications to polarizabilities in finite volume. These modifications are complicated by 
both 50(4) breaking and photon zero-mode interactions. We determine the finite volume 
modifications to pion Compton scattering. Here we argue, however, that these modifica- 
tions cannot be straightforwardly utilized to ascertain finite volume effects for background 
field calculations of polarizabilities in lattice QCD. Section HVl summarizes our work, and 
Appendix [B] collects the finite volume functions employed in the main text. 

II. PION COMPTON SCATTERING 
A. Compton Scattering Amplitude 

For Compton scattering in infinite volume, the amplitude for a real photon to scatter off 
a pion can be parametrized as 



expression, Q n is the charge of the pion in units of e > 0. In writing the physical amplitude, 




(1) 
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we have made use of Coulomb gauge in which the initial and final polarization vectors, e M 
and e'*, satisfy Eq = e'q = 0. The Compton amplitude appearing above, moreover, includes 
the one-particle reducible and irreducible pieces, as we have retained the Born terms. 

The frequency independent term proportional to Q\ reproduces the Thomson cross- 
section for low-energy scattering of charged particles when the amplitude squared is com- 
bined with appropriate phase-space factors. This term is exactly fixed by the total char ge o f 
the system in accordance with the Gell-Mann-Golberger-Low low energy theorems 0, 24 . 
The induced El-El interaction strength is the electric polarizability, while the induced 
Ml-Ml interaction strength (3m is the magnetic polarizability. In order to identify these as 
polarizabilities one must pull out a factor of twice the target mass from the Compton ampli- 
tude, as we have in Eq. ([I]). The electric and magnetic polarizabilities are the first structure 
dependent terms in the low-energy expansion of the Compton scattering amplitude. These 
polarizabilities can be determined from first principles using lattice QCD techniques. In 
order to make the connection between lattice data and real world QCD, extrapolations in 
quark mass and lattice volume are required. To perform these requisite extrapolations, we 
turn to the low-energy effective theory of QCD, x?T. 



B. PQxPT for Pion Compton Scattering 

In current lattice calculations, valence and sea quarks are often treated differently. In 
the rather extreme approximation known as quenched QCD, the sea quarks are completely 
absent. Less extreme is partially quenched QCD, where sea quarks are retained but have 
different masses than their valence counterparts. While both approximations are certainly 
contrary to nature, the latter retains QCD as a limit. Observables computed in partially 
quenched QCD can be connected to their real world values by utilizing partially quenched 
%PT (PQ^PT) to derive formulae for the required extrapolation in sea quark mass. Because 
%PT is contained as a limiting case of PQxPT, we focus our discussion on PQ^PT. Peculiar- 
ities of quenched %PT (Q%PT) will be noted where relevant and the general conventions ap 
pear in Appendix [X] For further details on QxPT and PQxPT, see Hi 0, M, 0, H M, E 



To determine pion observables, we imagine that the strange quark mass is fixed at the 
physical value so that no extrapolations are needed in the valence strange or sea strange 
quark masses. To this end, we consider a partially quenched theory of valence u and d 
quarks, paired with degenerate ghost quarks u and d, and two additional sea quarks j and 
I. The quark masses are given in a matrix 

rriQ = dmg(m u ,m d ,mj,mi,m u ,m d ) , (2) 

where the last two entries are ghost quark masses that are degenerate with their valence 
counterparts. For simplicity below, we work in the isospin limit of the valence and sea 
sectors, so that m u = m d and nij = mi. PQxPT describes the low-energy dynamics of 
partially quenched QCD and is written in terms of the mesons $ that appear in the coset 
field S as 1 

'2i& 



exp ( — ) . (3) 



In our conventions, the pion decay constant / « 132 MeV. 
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These mesons are the pseudo-Goldstone modes appearing from spontaneous chiral symmetry- 
breaking: SU(4\2)l <S> SU(A\2)r — > SU(A\2)v- The dynamics of these modes is described at 
leading-order by the PQxPT Lagrangian 

f 2 f 2 / \ 

£ = i-str (D^D^Y) + A^str (m^E + ^m Q j - m 2 $ 2 . (4) 

Here $o = str(<&) / \[2 is the flavor singlet field which has been included as a device to obtain 
the flavor neutral propagators in PQxPT. Expanding the Lagrangian to tree-level, one finds 
that mesons composed of a quark Qi and antiquark Q- have masses given by 

m Q z Qj = X \( m Qh + ( m Q)jj) • (5) 

The flavor singlet field additionally acquires a mass m 2 ,. In PQxPT (as well as in xPT), 
the strong U(1)a anomaly generates a mass for the flavor singlet field and hence mo can be 
taken on the order of the chiral symmetry breaking scale, mo ~ A x « 47r/. The flavor singlet 
field can thus be integrated out. Flavor neutral propagators in PQxPT, however, cannot be 
diagonalized into simple single pole forms [31[ . This fact notwithstanding, the results of our 
calculations will not require the explicit form of these flavor neutral propagators. 

In writing the above theory of mesons, we have added the effects of electromagnetism 
in the leading-order Lagrangian. The U(l) gauge field, A^, appears in the action of the 
covariant derivative, D^, namely 

= + ieA^ [Q, S] , (6) 

where Q is the quark electric charge matrix. To completely specify how electromagnetism 
is coupled, we must extend the quark charges to the partially quenched theory. The choice 

Q = diag (q u , q d , q j: qi, q u , Qd) , (7) 

with qj + qi ^ 0, is particularly useful because it retains sensitivity to all electromagnetic 
couplings in the theory as well as maintains the cancellation of disconnected operator inser- 
tions between the valence and ghost sectors [H, 33|. Other choices are possible but can be 



computationally cumbersome in actual lattice simulations. 



C. Pion Polarizabilities in Infinite Volume 

To determine the pion polarizabilities, we calculate the Compton scattering amplitude 
for the process 77T — > 77T using PQ^PT. Contributions to the amplitude are of three types: 
tree-level, wavefunction renormalization corrections, and one-loop contributions. The first 
contributions arise from tree-level graphs generated from local electromagnetic vertices in 
the effective theory. The relevant diagrams have been depicted in Figure [Tj, and are only 
non-vanishing for the charged pion. The first diagram represents the local coupling of two 
photons to the pion. This diagram arises from both the charge-squared operator contained in 
the leading-order Lagrangian, as well as from terms in the next-to-leading order Lagrangian. 



Specifically in the notation of [34( , the local two-photon, two-pion interactions are contained 



in the next-to-leading order terms 2 

C = iea 9 F flu str (QD^D U Y) + QD^D U E) + e 2 a 10 F 2 str (QEQTt) 



2 Although we use the SU(3) notation for these terms, final results depend on the scale-independent com- 
bination ag + aio, which has the same value in SU(2) as it does in 5/7(3). 
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FIG. 1: Tree-level contributions to the Compton scattering amplitude. The dashed lines represent 
mesons, while the wiggly lines represent photons. Vertices are generated from the leading and 
next-to-leading order Lagrangian. 




FIG. 2: Wavefunction renormalization in PQxPT. Diagram elements are the sames as in Figure [TJ 
and the cross denotes the partially quenched hairpin. The vertex is generated by the leading-order 
Lagrangian. 



where F^ u = d^A v — d v A^ is the electromagnetic field-strength tensor. In PQxPT, the 
low-energy constants a>g and a\o have the same numerical values as in xPT, which can be 
demonstrated by matching. The remaining two diagrams in Figure [1] are Born terms that 
do not contribute to the one-pion irreducible Compton amplitude. 

The next contributions are those that arise from the pion wavefunction renormalization. 
The leading self-energy diagrams are depicted in Figure [21 The leading-order diagrams 
involving the photon coupling to the pion charge must be multiplied by the wavefunction 
renormalization to obtain contributions to the Compton amplitude at next-to-leading order. 
Thus we require only the wavefunction renormalization of the the charged pion. Due to 
fortuitous cancellation in both PQxPT and QxPT, the hairpin diagram, which arises from 
the double pole structure of the flavor-neutral propagator, vanishes. 

The remaining contributions to the Compton amplitude arise from one-loop diagrams. In 
Figure [31 we display the diagrams for the one-pion irreducible scattering amplitude. Contri- 
butions from such diagrams lead to chiral corrections to the electromagnetic polarizabilities. 
For the charged pion, there are additional one-pion reducible pieces in PQxPT. These dia- 
grams are displayed in Figure HI The effects of such diagrams in infinite volume, however, 
are to renormalize the mass of the intermediate state pion, and to provide the necessary 
cancellations which preserve the charge interaction of the leading Born terms. The latter 
cancellations were first worked out explicitly for the case of the pion charge radius in PQ^PT 



m 



361 ] . Assembling the results of Figures [TH3 we can extract the pion polarizabilities 
using Eq. ([1]) by utilizing Coulomb gauge in the center-of-momentum frame. 

At one-loop order, it is well known that (Xe+Pm = for both charged and neutral pions [!, 



371 . |38| . We find this remains true to one-loop order in PQ^PT, as well as Q%PT. This is 
expected because extending the flavor algebra from SU(2) to graded Lie algebras cannot 
alter the helicity structure of the Compton amplitude. As for the orthogonal combination 
of polarizabilities, a^- (3m, we arrive at 

*° on 2a fsQl fQ S 

n«± _ lGa f8 Ql{a Q + ai Q ) , , 

with Q n = q u — q d . These results are the same in xPT, PQxPT, and QxPT, with the 
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FIG. 3: One-loop contributions to the Compton scattering amplitude in PQxPT. Vertices shown 
are generated from the leading-order Lagrangian, and diagrams depicted are all one-pion irre- 
ducible. 



exception that in QxPT the low-energy constants ag + ctio, and / have different numerical 
values. Furthermore our x?T result agrees with the literature, see (being careful to note 
/ = y/2f w ). I n deriving the above result, we remark that the delicate cancellation between 
pion loops in the zero frequency limit present in %PT remains in PQxPT, and QxPT. 
This cancellation is required by the infinite volume gauge invariance of the amplitude and 
reflects that the longest-range coupling to the pion is only to the total charge. In this 
way the Thomson scattering cross section is produced in these three theories when the zero 
frequency limit is taken. 

In each theory, there are no local electromagnetic interaction terms for the neutral pion in 
the next-to-leading order Lagrangian. Thus there can be no divergences in the polarizabilities 
of the neutral pion, as we found explicitly at one-loop. While chiral logarithms are absent 
for the neutral pion, there are finite terms from the loop graphs. In fact, the entire pion 
cloud contribution to the neutral pion polarizabilities manages to survive quenching. This is 
rather surprising, but can be understood by considering the quark-line topologies generated 
at one-loop order. 

The five topologies arising from the four-pion vertex generated from Eq. (j3J) are depicted 
in Figure The topologies in the second row are only possible for flavor neutral external 
states, such as the neutral pion. Let us investigate which topologies can make non-vanishing 
contributions to the neutral pion polarizabilities. Diagram A contains a sea quark loop and 
thus associated contributions are proportional to 



which sums the charge-squared couplings from all possible valence-sea loop mesons. In the 
isospin limit of S77(4|2) all such mesons are degenerate with mass rrij U . The net contribution 
from topology A, however, vanishes because contributions from the four-meson vertex with 
two derivatives exactly cancel contributions from the four-meson vertex with quark mass in- 



AQ 2 = (q u - qj f + (q u - qif + (q d - qj) 2 + (q d - q t ) 2 , 



(11) 
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FIG. 4: One-pion reducible contributions to the Compton scattering amplitude at one-loop order 
in PQxPT. 




FIG. 5: Quark-line topologies generated at one-loop order. Diagrams in the second row contribute 
only when the external states are flavor neutral. 

sertion. Terms from all of the meson loop diagrams in Figure [3] are required for this delicate 
cancellation. As a result, the characteristic factor of AQ 2 is absent from Eq. ([9]). Next, each 
of the quark line topologies B, C, and D, require flavor disconnected contributions from 
flavor-neutral meson propagators. As flavor neutral mesons are also electrically neutral, 
coupling to the photon eliminates such contributions. Indeed looking at Figure El there is 
only one possible diagram with a hairpin vertex. Direct evaluation shows that this contri- 
bution vanishes, ruling out the B, C, and D topologies. Therefore the loop contributions to 



S 



neutral pion polarizabilities in Eq. (Q stem entirely from topology E. As this topology is 
quark-line connected, the contribution has the same form regardless of quenching. 

Let us examine topology E closer by writing the pion field in terms of quark basis mesons, 
7] u and ri d , 

7r° = ^=( Vn - Vd ). (12) 

The diagonal contractions, r] u -r] u and r] d -T] d , for topology E result in an electrically neutral 
loop meson. Only the first diagram of Figure [3] could yield the diagonal contractions of 
topology E. Close inspection of the Lagrangian shows, however, that the four-meson, two- 
photon vertex with four electrically neutral mesons is identically zero. Thus neutral pion 
polarizabilities stem entirely from topology E's non-diagonal flavor contractions: Tju-rjd, and 
Vd-Vu- To one-loop order, the neutral pion polarizabilities arise entirely from annihilation 
contractions of lattice QCD correlation functions. 

Returning to Eq. ( fTOl) for the charged pion polarizabilities, the quark-line picture helps 
to show why Eq. ([TO]) has no loop contributions. In Figure the quark-line topologies in 
the second row are no longer relevant because the external states are charged. Furthermore 
topologies B and C require hairpins, but the hairpin graphs in Figures [2] and [3] vanish. Loop 
contributions to charged pion polarizabilities can only arise from topology A. Cancellation 
of divergent loop contributions from this topology must occur in xPT, PQ^PT, and QxPT 
because the combination a$ + ctio is renormalization scale independent. This independence 
disallows chiral logarithms from loop contributions in %PT, and QxPT. While one can 
imagine scale invariant combinations of the form log(m?/ rrif u ) , say, away from the isopsin 
limit of S77(4|2), charge- squared couplings do not allow for loop contributions to alternate 
in sign, see Eq. (ITT]) . Thus such logarithms are absent. While logarithms are not allowed, 
finite contributions can be present. As with the neutral pions for topology A, however, the 
contributions from the four-meson vertex with two derivatives exactly cancel contributions 
from the four-meson vertex with quark mass insertion. The characteristic factor of AQ 2 
is consequently absent from from the charged pion polarizabilities, Eq. (TT0|) . Thus the 
accidental cancellation of finite terms in xPT also occurs in PQ^PT, and QxPT. 

As a final comment on the inifinite volume results in Eqs. and (1101) . the only pion 
mass dependence in both charged and neutral pion polarizabilities arises from the target 
mass m^. The target mass depends on the valence quark mass. Both of these statements 
hold only to one-loop order in the chiral expansion. 



III. COMPTON TENSOR IN FINITE VOLUME 

In finite volume, the pion is already deformed. Thus its ability to polarize in an applied 
electric or magnetic field will differ from that in infinite volume. As finite volume modifica- 
tions to hadron properties are long distance in nature, they can be quite generally addressed 
using xPT. Lattice simulations are usually carried out in a hypercubic box of volume L 3 x (3, 
where L is the length of the spatial direction, and (3 is the length of the Euclidean time di- 
rection. We consider (3 ^> L so that there is no effect from the finite temporal extent of 
the lattice. With periodic boundary conditions on the quark fields in each of the spatial 
directions, the momentum modes on the lattice are p = 2irn/L, with n a triplet of integers. 
The ordinary power counting for ;\PT 

\p\ < m n < A x (13) 
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can be applied in a box of finite size provided 2fL ^> 1 and m^L > 2n. These two conditions 
constitute what is called the p-regime of chiral perturbation theory. The first condition 
is required in order to use chiral perturbation theory at all, while the second condition 



maintains that pionic zero modes remain perturbativ e |39l . 401. As the power counting in 



this regime remains the same as in infinite volume [4l|], the same diagrams depicted in 
Figures [TH contribute to the pion Compton tensor. It is straightforward to perform the 
loop calculations in a finite box, we merely replace integrals over virtual four-momenta by 
integrals over energy and sums over spatial momentum modes permitted by periodicity. 

Consider an observable X calculated in both finite and infinite volume. Let X(L) denote 
the value of the observable in finite volume, and X(oo) denote its value in infinite volume. 
The finite and infinite volume theories share exactly the same ultraviolet divergences, thus 
the volume effect can be determined from matching the two theories in the infrared, 

X(L) =X(oo) + AX(L). (14) 

The volume effect is given by the matching term AX(L) which is ultraviolet finite. A salient 
feature of such matching is that it allows us to retain our infinite volume regularization 
scheme and values of low-energy constants. 

Calculating the finite volume matching for the Compton scattering amplitude, while 
straightforward, is quite involved in practice. We first remark that the decomposition of the 
Compton tensor in Eq. ([T]) is no longer valid. That decomposition makes use of the center- 
of-momentum frame. Finite volume results on a torus have only an S4 cubic subgroup of 
the infinite volume S"0(4) invariance. Results will thus be frame dependent, and hence, to 



be general, we must not make recourse to a particular frame. Furthermore, as shown in [42 
there are more gauge invariant structures allowed on a torus. Thus more terms than shown 
in Eq. ([1]) are allowed at second order in the field strength. 
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A. Neutral Pion 



Carrying out the finite volume matching on the neutral pion Compton amplitude without 
recourse to a particular frame or gauge, we find 



gt*" / dxl 3/2 (xr, m% - x(l - x)r z ) 



1 



+5 ^S v0 / _ J dx J d yl 3/2 (xk + yk', ml - xyr 2 ) 



1 



o Jo 



dx / dy[(2x - l)ui + 2yuj'][2xw + (2y - l)u']I 5 / 2 {xk + yk' 1 m 2 tj) - xyr 2 ) 



jA^S"! [ dx I dy[(2x~l)uj + 2yuj'}\k' j I 5/2 (xk + yk',ml-xyr 2 ) 
4 Jo Jo 

+ 2I J 5/2 (xk + yk', ml - xyr 2 ) j 

+ -S^S v0 [ dx [ dy[(2y- 1)lj' + 2xoj]\k i I h/2 (xk + yk\nA -xyr 2 ) 
4 Jo Jo 

+ 2I* /2 (xk + yk',ml-xyr 2 )} 

i r 1 r l ~ x 

_ -Si»S" 3 I dx dy[4f 5 :l /2 (xk + yk',ml-xyr 2 ) + 2k l P 5/2 (xk + yk',ml-xyr 2 ) 

4 Jo Jo 

+ 2k' j ll /2 {xk + yk', ml - xyr 2 ) + Wk' 3 ' I s/2 (xk + yk' , m| - xyr 2 )] 



(15) 



Above we have employed r M = (k — k')^ for the momentum transfer, with k^ = (ou, uik), and 
k'^ = (u',u'k') for the photon momenta. To derive the above result, we have employed the 
reality of initial and final state photons, and taken their spatial momenta to be quantized. 
The finite volume functions 1/3(0, m 2 ), I%(0, m 2 ), and 1^(0, m 2 ) are defined in Appendix iBl 

The coefficient for contributing loop mesons CJ° is given by 



Cf = 3Q 2 n (2ml - r 2 ) <W - -AQ 2 r%, ju . 



(16) 



While we have only given the PQxPT coefficients, the xPT, and QxPT results can be 
trivially deduced from Eq. ffTB"]) . The latter is possible because there are no hairpin contri- 
butions. At zero frequency, the finite volume Compton amplitude for the neutral pion is 
non-vanishing. This is because the Thomson cross-section is not protected from renormal- 
ization in finite volume 14211 . 



B. Charged Pion 

The charged pion Compton amplitude at finite volume is even more involved than the 
neutral result as we must determine both reducible and irreducible contributions. To one- 
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loop order, the result is 

3e 2 AQ 2 2 



AT^(L) 



2/2 

+<j"V° f dx [ 
Jo Jo 



1 f 1 

q9* v J o dxl 3/2 (xr,'. 



dy 



~-I 3 / 2 {xk + yk', m 2 u - xyr 2 ) 



-[(2s - l)w + 2yu'][2xw + (2y - l)u>']I 5/2 {xk + yk' , m 2 „ - xyr 2 ) 



+ -5^S vj I dx 



o Jo 



l-x 



dy[(2x - l)oj + 2yJ\ 



k' J I 5 / 2 (xk + yk 



' m 2 



xyr 2 ) + 2I 3 5/2 (xk + yk' , m 2 u - xyr 2 ) 



+ -5^5"° 
4 



l-x 



dy[(2y - l)u/ + 2xw] 



dx 
Jo 

k l I 5/2 (xk + yk', m 2 u - xyr 2 ) + 2P 5/2 (xk + yk', m 2 u - xyr 2 ) 



-Isnipi J dx J dy [4ll j /2 (xk + yk', m 2 u - xyr 2 ) + 2k% /2 (xk + yk' , m 2 u - xyr 2 ) 
+ 2k' j P 5/2 (xk + yk', m 2 u - xyr 2 ) + tfk' 3 I b/2 (xk + yk' , m 2 u - xyr 2 ) \ 



2e 2 Qj 
f 2 

-e 2 Ql 
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(2P + fc)"(P + P' + k) v (2P - k') v (P + P' - k') 



(P + k) 2 - ml 



(P - k') 2 - ml 



h/2{m 2 ju ) 



I»{P,P + k)(P + P' + k) v {2P + kYI v {P + k,P') 



(P + k) 2 - ml 



(P + k) 2 -ml 



i»(p - k', p'){2P - k'f (p + p 1 - k'yr{p, p - k') 



(P-k') 2 -ml 



(P - k') 2 - ml 



+ 



e 2 Ql 



S^S™ I dx 
o 



§ vj f 

Jo 

f 

Jo 



dx 
dx 



-6^6 VJ / dx 
Jo 



2I 1/2 (xk,m ju ) + 2I 1/2 (xk ,m ju ) 
x(2x - 1) (uj 2 I 3/2 {xk, m 2 u ) + uj' 2 I 3/2 {xk' , m 2 J) 
xu'k'n 3/2 {xk' ,m 2 u ) + 2xu'P 3/2 (xk',m 2 u ) + (2x - l)uP 3/2 {xk,m] 

xuk l I 3j2 (xk, m 2 u ) + 2xivP 3/2 (xk, m 2 u ) + (2x - l)JP 3/2 {xk', m 2 u ) 
2P 3 3 l2 {xk, m-J + 2Il 3 /2 (xk', to 2 J 



+ k l P 3/2 (xk, to 2 J + k ,j P 3/2 {xk', m 2 u ) 



(17) 



In the above result, the initial (final) pion momentum has been denoted by P (P 1 )- We have 
employed an abbreviation for the finite volume pion-photon vertex function, I fl (P2, Pi), 
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which arises from the one-pion reducible diagrams and is given by 



nP 2 ,Pi) = 72^°j( P 2 + A)°[y o dxI 1/2 (xA 7 m 2 u -x(l-x)A 2 ) 



- -(A°) J dxx(2x - l)/ 3/2 (xA, mj u - x(l - x)A 2 ) 
^A° J dx(l - 2x)(P 2 + A) ■ I 3/2 (xA, m% - x(l - x)A 2 



5M{{P 2 + P 1 ) 1 / dxr 3 j /2 (xA,m 2 ju -x(l-x)A 2 ) 



\tf I dxP 3/2 (xA,m 2 u -x(l~x)A 2 ^ 
1 



A°(P 2 + P 1 )° 
1 



dxx 



I J 3/2 (xA,m 2 u ~x(l~x)A 2 



2 A^ y dx/ 3/2 (xA, m 2 ju - x(l - a;)A J 



(18) 



with A M = (P 2 — Pi) 11 - At zero frequency, we recover the results of 42|. Specifically from 
the one-pion reducible terms, we see that the current is renormalized. This is possible at 
finite volume because of gauge invariant zero-mode interactions. 



C. Discussion of Finite Volume Results 

With Eqs. ( fT5l) and (TTTI) . we have deduced the finite volume modification to the pion 
Compton scattering tensor. These results show explicitly broken 50(4) invariance as well 
as additional structures not anticipated by infinite volume gauge invariance. The finite vol- 
ume modifications can be directly utilized if two-current, two-pion correlation functions are 
calculated on the lattice. One merely removes the finite volume effects determined above 
to isolate the infinite volume physics. Such lattice calculations of the Compton tensor are, 
however, prohibitively expensive time wise, and will not be performed in the foreseeable 
future. A practical alternative to these calculations is provided by the background field 
method. In this approach, a classical electromagnetic field is gauged into the QCD action. 3 
One then studies the external field dependence of correlation functions to deduce electro- 
magnetic observables. For example, at infinite volume the energy of a neutral pion in a weak 
external electric field is 



E n (p = 0) = m n - l -alE 2 + 0{E A ) 



(19) 



Thus by measuring the quadratic energy shift in the external field strength \E\, one can 
deduce the electric polarizability. A practical question is then how to deduce volume cor- 
rections to polarizabilities determined from background field methods. Given the relation 
of the infinite volume Compton tensor to the polarizabilities, one might suspect that the 



3 Implementing this method currently suffers the need to quench effects of the background field. In principle, 
there is no impediment to coupling a suitably weak background field to sea quarks other than time cost. 
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finite volume Compton tensor in Eqs. (|T5|) and (jTTj) contains the finite volume corrections 
to the polarizabilities. We argue that the finite volume Compton tensor has no relevance to 
volume effects in background field methods. At finite volume, there is no longer a discernible 
relation between polarizabilities and the Compton tensor. 

An analysis of finite volume effects for nucleon polarizabilities for background field meth- 



ods derived from the Compton tensor, however, was presented in [21[. That analysis em- 
ployed the Breit frame decomposition of the nucleon Compton tensor in Coulomb gauge. 
Such results surely cannot be utilized for background field calculations because such cal- 
culations are typically done in the rest frame. The finite volume modifications derived, 
moreover, are polluted by subtle effects from the gauge field due to the nature of gauge 
invariance on a torus. These effects have nothing to do with polarizabilities. For example, 
in the center-of-momentum frame, where k° = k'° = lu, we may encounter a term in the 
amplitude of the form 

M — ...+ \^ 2 a{L) e'* ■£ + ..., (20) 

and be tempted to conclude that a(L) is a finite volume correction to the electric polariz- 
ability. In a general frame, however, this term could stem from any combination of u 2 , a/ 2 , 
and uouj' structures. In infinite volume only the last term is allowed by gauge invariance, 
specifically by an operator oc E 2 with a coefficient proportional to the electric polarizability. 
In finite volume, however, the additional structures u 2 and a/ 2 are allowed. They stem from 
single-particle effective theory operators of the form 

£= l -a(L)W^-^tr(Q 2 <t> 2 ), (21) 

for example, where is the negative parity part of the zero-mode Wilson line Wi, given 

by 

Wi = VoWiVl (22) 

with the Wilson line W; as 

W t = exp r^jdxiA^j , (23) 

and Vq as the zero-mode projection operator. The operator in Eq. ( 12T1) respects C, P, and 
T, as well as the cubic symmetry of the torus. Furthermore it is gauge invariant because the 
zero mode has a periodicity constraint under gauge transformations, see j42|. From Eq. (1201) . 
we cannot deduce that a(L) is a finite volume correction to the electric polarizability. In 
general, one must work in an arbitrary frame to disentangle the zero-mode electric coupling 
in Eq. (12 ip from the electric polarizability. An analogous situation exists for magnetic 
interactions, because the operator, 

V • (W { - } x B) tr(Q 2 $ 2 ), (24) 

for example, is allowed by symmetries. 

The frame and gauge dependence notwithstanding, finite volume modifications to polar- 
izabilities were determined in (2l| from Taylor series expanding the Compton amplitude in 
photon frequency. That procedure is also invalid as we now demonstrate. For simplicity, 
consider the following finite volume difference function, /^(A^m 2 ), where k is an exter- 
nal photon momentum. To determine finite volume corrections to polarizabilities stemming 
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from this term, we perform a Taylor series expansion in the external momentum and arrive 
at 

h /2 (k,m 2 ) = / 1/2 (0,m 2 ) - ifc 2 m 2 / 5/2 (0,m 2 ) + 0(k A ). (25) 

If we were interested in determining a hypothetical polarizability X entering the amplitude 
in the form 

M = ... + ^k 2 X + (26) 

then we would be tempted to conclude that the finite volume effect AX(L) is given by 

AX(L) = -m 2 / 5 / 2 (0,m 2 ). (27) 

Because the external momentum is itself quantized, instead of Eq. (l25l) we actually have the 
exact relation 

J 1/2 (fc,m 2 ) = / 1/2 (0,m 2 ). (28) 

This follows trivially from re-indexing the summation over loop momentum modes, or from 
the periodicity of the elliptic-theta function, see Appendix [B] Hence the volume effect for 
our example is actually AX(L) = 0. The reason for this discrepancy is a poorly convergent 
series expansion. 4 Naively the expansion is in k 2 = 4ir 2 n 2 /L 2 , and thus for large enough 
box size the finite volume effect should be well approximated by the first few terms in the 
Taylor series. This is not the case. Because higher-order terms have more derivatives, these 
contributions effectively have more propagators and hence more sensitivity to the infrared. 
While we would expect the second term in Eq. fl25]) to be 1/L 2 suppressed relative to the 
first term, the asymptotics show that the volume effect is L 2 enhanced 

Km m 2 J 5/2 (0,m 2 )// 1/2 (0,m 2 ) = \h 2 . (29) 

L — >oo O 

The series expansion continues in this fashion: all terms are order one. We can see the same 
effect more directly by expressing the finite volume difference in terms of the elliptic-theta 
function, namely 

h/2(k, m 2 ) = / d\e- m2L2 ^ x # 3 (™, e- A )tf 3 (0, e" A ) 2 - 1 , (30) 

for the choice k = (27m/ L, 0, 0). A series expansion in k is thus effectively the same as 
expanding in irn. 

Returning to Eqs. (|15p and (|17|) . we must ascertain whether we can make sense of a series 
expansion in frequency for the Compton tensor. Terms of the form 

r 2 / dxl 3 / 2 (xr, m 2 — x(l — x)r 2 ), (31) 
Jo 

for example, can be plausibly expanded to second order because this requires only evaluation 
of the finite volume function at r = 0. This was the logic employed in [431 to deduce finite 



4 Another difference between Eqs. (|25|) and (|28|) is that the order of summation and differentiation has been 
interchanged. One can easily show, however, that the summation over modes is uniformly convergent by 
using the Weierstrass M-test. 
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volume corrections to the nucleoli magnetic moment. As r is not continuous, however, one 
cannot deduce the small momentum behavior of this term from evaluation at r — 0, nor can 
one deduce the small momentum behavior from Taylor series expanding other terms like 



in the Compton tensor, for example. Series expanding in rL = 2irn is nonsense no matter 
the size of the box length L. 5 

At finite volume, one must treat the terms in the amplitude as form factors in uoL. 
Thus for electromagnetic form factors at finite volume, for example, volume corrections to 
electromagnetic moments cannot be deduced. Similarly we are unable to use our results for 
the finite volume modification of the Compton amplitude to deduce corrections to the pion 
polarizabilities. At second order in the field strength there are a myriad of new terms allowed 
by the less restrictive symmetries on a torus: cubic invariance and periodic zero-mode gauge 
invariance. Furthermore a small frequency expansion at finite volume does not make sense 
for quantized momenta. Said another way, periodic gauge potentials on a torus do not lead 
to electromagnetic multipole expansions. 

IV. SUMMARY 

Above we have investigated chiral and volume corrections to pion Compton scattering us- 
ing xPT, PQxPT, and QxPT. In infinte volume, straightforward calculation of the Compton 
amplitude allows us to determine charged and neutral pion polarizabilities in these theories. 
Due to fortuitous cancellation there is no dependence on the sea quark masses, or sea quark 
charges at one-loop order in the chiral expansion. The Compton tensor itself does not have 
any quark mass dependence at this order. Consequently the quark mass dependence of the 
derived polarizabilities stems from a kinematical prefactor of the inverse target mass. As 
this valence pion mass is relatively inexpensive to dial, the chiral singularity should be dis- 
cernible from lattice data at light quark masses. Thus as the chiral regime is approached, 
one can use the lattice as a diagnostic tool to study the chiral behavior of pion polarizabil- 
ities. This can be done most easily for the charged pion. Whereas for the neutral pion, 
we demonstrated that the polarizabilities at one-loop order stem entirely from annihilation 
contractions which are notoriously difficult to calculate on the lattice. 

When accounting for finite volume effects, however, the situation becomes more compli- 
cated. Breaking of 5*0(4) invariance and the nature of gauge invariance on a torus lead 
to considerably complicated structure for the Compton tensor. Sea quark charge and mass 
dependence enter in the one-loop finite volume effects. One cannot unambiguously deter- 
mine the volume effects for the polarizabilities from the Compton tensor because the Taylor 
series expansion in quantized momentum is poorly convergent. What was in infinte volume 
a series expansion in u/m n <C 1 that lead to the polarizabilities, now is accompanied by 
an ill-defined expansion in uL ~ 1 at finite volume. This means that even at low energies, 
the finite volume Compton amplitude is a form factor in uL. Consequently connection of 
our finite volume results to background field lattice calculations is not possible. Similarly 



5 There is a putative improvement in the convergence of the last term due to the integral over the Feynman 
parameter. The L scaling of terms in the expansion, however, is unchanged. 




(32) 
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finite volume corrections to electromagnetic moments cannot be deduced from momentum 
expanding finite volume form factors. Further investigation is required to determine volume 
corrections relevant for observables determined with background field methods. 
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APPENDIX A: QUENCHED X PT 

Here we give the relevant details needed in the calculation of quenched pion polariz- 
abilities. In quenched QCD, contributions from sea quarks are completely neglected^ In a 
quenched theory of two flavors u and d, we additionally have two ghost quarks u and d. The 
mass matrix is now 



where the final two entries are the masses of the ghost quarks. These equal mass ghost 
quarks are necessitated so that path integral determinants for the valence quarks are exactly 
canceled by those from the ghosts. The symmetry breaking pattern in Qx?T schematically 
takes the form C/(2|2)l<8)[/(2|2).r — > U(2\2)y because there is no axial anomaly in quenched 
QCD. The coset field £ is hence a U{2\2) matrix and the singlet component cannot be 
integrated out. The dynamics of the pseudo-Goldstone modes is described at leading-order 
by the Q%PT Lagrangian 



While propagators for flavor-neutral mesons have double poles, these are not encountered 
explicitly in expressions for the pion polarizabilities at next-to-leading order. Quenched 
observables are in general unrelated to their unquenched counterparts, for example, the 
constants a$ and mo have no analogs in xPT, moreover, arbitrary polynomial functions of $q 
can multiply any term in the Lagrangian and the low-energy constants in the quenched chiral 
Lagrangian above result from treating these polynomial terms in mean-field approximation. 
For the quenched electric charge matrix of the quarks, we must have 



for which the condition strQ = is unavoidable. In general there are fewer local electromag- 
netic terms in Q%PT as compared to %PT. At next-to-leading order, however, both the ag 
and aio terms remain. We must keep in mind that the numerical values of these coefficients 
are unrelated to their values in x?T. Calculation of the pion polarizabilities then proceeds 
analogously to the partially quenched case. Results for the quenched polarizabilities have 
been given for infinite and finite volume in the main text. 



m Q = diag (m u , m d , m u , m d ) , 



(Al) 




Q = diag (q u , q d , q u , Qd) , 



(A3) 
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APPENDIX B: FINITE VOLUME FUNCTIONS 



Above we have determined the finite volume modification to the Compton scattering 
tensor. In this Appendix, we give explicit formulae for the finite volume functions used to 



express finite volume differences. We use similar notation for these functions as |44j . |45 
where further discussion can be found. 

In evaluating a Feynman diagram in finite volume, the loop integral is converted into a 
sum over the allowed Fourier modes in a periodic box. The difference of this sum and the 
infinite volume result is the finite volume effect. As is customary, we treat the length of 
the time direction as infinite. All finite volume differences with momentum insertion can be 

defined by 



cast in terms of the function 1^ 3 (0, m 



1 



qii . . . q l i 



L 3 ^ [(q + 0) 2 + m 



21/3 



dq 



q n ■ ■ ■ q l i 



(2vr) 3 \(q + Of + m 2 }^ 



(Bl) 



where the sum on n is over triplets of integers, and the loop momentum modes are quantized 
as q = 2irn/L in a periodic box. While a general expression for the exponentially convergent 

form of %J (6, m 2 ) exists, it is easiest merely to cite the required cases for our work. These 
are 



%(6, m 2 ) 



(L 2 /4) 



(3-3/2 roo 



(47T)3/2r(/3) J 
1 



rfAA 1 / 2 - /3 e" m2i2 / 4A 



d 



3=1 
2\ 



2{P - 1) d&- 
1 



-Ip{0,m')-6^Ip{0,m l ) 
d 2 



A((3 -2)((3 -l)de^d6^ 



1(3-2(0, m 2 ) 



2((3-l 



- e h r f3 2 (0,m 2 ) - e i2 f^{6,m 2 ) - e ii e i2 i l3 {6, m 2 ) 

where $3(2, q) is a Jacobi elliptic-theta function of the third kind, see, e.g. [46]. 
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